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Abstract The steady equilibrium conditions for a mixed gas of neutrons, protons, elec- 
trons, positrons and radiation field (abbreviated as npe ± gas) with/without external 
neutrino flux are investigated, and a general chemical potential equilibrium equation 
[i n — pL p + Cp e is obtained to describe the steady equilibrium at high temperatures 
(T > 10 9 K). An analytic fitting formula of coefficient C is presented for the sake of sim- 
plicity as the neutrino and antineutrino are transparent. It is a simple method to estimate 
the electron fraction for the steady equilibrium npe ± gas that using the corresponding 
equilibrium condition. As an example, we apply this method to the GRB accretion disk 
and approve the composition in the inner region is approximate equilibrium as the accre- 
tion rate is low. For the case with external neutrino flux, we calculate the initial electron 
fraction of neutrino-driven wind from proto-neutron star model M15-ll-rl. The results 
show that the improved equilibrium condition makes the electron fraction decrease sig- 
nificantly than the case fi n = fj, p + p e when the time is less than 5 seconds post bounce, 
which may be useful for the r-process nucleosynthesis. 

Key words: nuclear reactions, nucleosynthesis, weak-interaction, GRB accretion disk, 
neutrino-driven wind 



1 INTRODUCTION 

It is a classic and simple approximation for the practical application at many astrophysical sites that 
matter compositions can be considered as a mixture of the neutrons, protons and electrons, i.e. so called 
npe~ system. If the temperature is very high(T > 10 9 K), lots of photons, positrons, even neutrinos 
and antineutrinos will appear in the system, i.e., the system becomes a mixture of electrons, positrons, 
nucleons and radiation field (we abbreviated it as npe* gas). Many astrophysical sites can be regarded 
as the n pe^ gas, such as (i) th e hot fireball jetted from a successful central engine of Gamma Ray Burst 
(GRB) dPruet & Dalall |2002|) . (ii ) the matter after the c ore-collapse supernova shock due to the photo- 
disintegration of the iron nuclei dMarek & Jankal 120091). (iii)th e neutrino-driven wind comes from the 
prot o-neutron star(PNS) as T > 10 9 K qM artmez-Pinedo, 2008), (iv) the outer core of t he young neutron 
star dYakovlev et all 120081: iBaldo & Ducoinll2009h . ("v)tne accreting disk of the GRBs dLiu et al.[ [20071; 
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Janiuk & Yuanll2010h and (vi) the early universe before the decoupling of neutrinos dDuttaet all 12004 
HarwitLl2006lT [n a word, npe and npe^ gas is applied widely up to the present. Steady equilibrium 
state of npe~ or npe ± gas is an important stage for many cases. Many authors have addressed this issue 
for several decades. A typic al disposal to the steady sate equilibrium of npe~ system is concluded by 
IShapiro & Teukolskvl ( ll983l) . They gave an important result that /i n = fj, p + fi e for a steady equilibrium 
npe" system, where p s are the chemical potentials for neutron,proton and electron respectively. This 
result have been accepted by most authors. But in fact, Shapiro et al. only considered the electron cap- 
ture and its reverse interaction at 'low tempera ture', they ig nored the appearance of positrons when the 
temperature of system is high enough. Recently ! Yuan! d2005l) argued that lots of positions can exit at high 
temperature, which leads to the great increase of the position capture rate. The positron capture affects 
the condition of steady equilibrium significantly. If the neutrinos can escape freely from the system with 
plenty of e ± pairs, the equilibrium condition should be /i n = p, p + 2p e instead. However, for a more 
general condition when the temperature is moderate, the equilibrium condition have not be researched. 
Liu et al. have ever taken a method in which they assume the co efficient o f p e va ries exponentially from 
Pn = Mp + Me to Mn = Pp + in the accreting disk of GRBs dLiu et al .1120071) . but it is not a rigorous 
method. Therefore a detailed and reliable database or fitting function to describe steady equilibrium of 
the npe^ 1 gas at any temperature is necessary. Furthermore, the above discussions are limited to the iso- 
lated system, ignoring the external neutrino flux. In this paper we investigate the chemical equilibrium 
condition for npe gas at any temperature from 10 9 to 10 11 K, and give a concrete application to the 
GRB accretion disk. We also calculate the initial electron fraction of the neutrino-driven wind in PNS, 
in which the external strong neutrino flux can not be ignored. This paper is organized as follows. In 
section II, we present the equilibrium conditions as neutrino is transparent or opaque for an isolated 
system. Section III contains a detailed discussion to the initial electron fraction of neutrino-driven wind 
for PNS model M15-ll-rl. Finally we analyze the results and make our conclusions. 



2 EQUILIBRIUM CONDITION OF NPE ± GAS WITHOUT THE EXTERNAL NEUTRINO 
FLUX 

For a mixed gas of npe ± and radiation field at different physical conditions, we divide them into 
two cases: neutrino transparence and opacity. To guarantee the self-consistence, we give a simple 
estimate for the opaque critical density of the npe ± gas. The mean free path of neutrino is l v = 
n<r sac +n a°- b3 ' wnere 71 ar, d n n are the number density of baryon and neutron respectively, n = pNa, 
n n = p(l — Y e )Nji, p is the mass density, Y e is the electron fraction, Na is the Avogadro's con- 
stant, af, ac and af s are the scatter cross section with b aryons and absorption section by the neutrons. 
a sac w (-Sj T ) 2 lQ- 44 (|Kippenhanhn & Weigeri Il990l) . where E v is the energy of neutrino, m e c 2 is 

the mass energy of electron, c is the light velocity. a® bs « -§iE e p e as -4-E^ dOian & Wooslevl 1 19961: 
lLai&Oianl[T998l) . where A = nG 2 F cos 2 c {C v + 3C\), G F = 1.436 x 10" 49 erg cm 3 is the Fermi 
weak interaction constant, cos 2 £? c = 0.95 refers to Cabbibo angle. Cy = 1, Ca = 1.26, E e andp e are 
the energy and momentum for electron, respectively. Due to the energy conversation of nuclear reaction, 
E e = E v + Q, Q = (m n — m p )c 2 = 1.29MeV, m n and m p are the mass of neutron and proton. At high 
density, the electrons are strong degenerate and relativistic, so E e w Ep = [(37r 2 Ag7i e ) 2 / 3 + l] 1 / 2 ^ 
unite of m e c 2 ), A e = is the reduced electron Compton wavelength. Substituting pY e NA for n e , 

E e w (3ir 2 \ 3 e pY e N A ) 1/3 - Therefore, the mean free path of neutrino is 



" pN A [(.3n 2 ~XlpY e N A - Q) 2 ' 3 x lO" 44 ] + [^(37r 2 ~\lpY e N A ) 2 / 3 }p(l - Y e )N A ' 

If we assume l v = 10km is criterion of neutrino opacity, p v cri = 5.58, 4.50, 4.10, 3.96, 3.96 x 10 10 g 
cm -3 for Y e = 0.1, 0.2, 0.3, 0.4, 0.5. Rigorously speaking, here we overestimate the absorption section 
because we ignore a block factor (1 — / e ), so p v cri is the minimum for critical density. As p < p v cri , 
neutrino is transparent, or it is opaque. By similar method, one only needs to replace v, E e = E v + Q, 
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and n„ to v, E e = E p — Q, and n p respectively for mean free path of antineutrino. The critical density 
for antineutrino p" crl = 1.43, 0.86, 0.62, 0.48, 0.40 x 10 n g cm" 3 for Y e = 0.1, 0.2, 0.3, 0.4, 0.5. 

Another preciser way to judge the transparency of neutrino is defining a parameter: neu- 
trino optical depth t, which is c l osely relative to object's composition and structure, r = 
< n e ff > (ir dArcones et all H008), where r is neutrino transport distance, < n e jf >= 

V< Kabs > (< Habs > + < K sac >), k q (, s and K sac are the absorption opacity and scatter opacity, 
Ksac — na sac , K a b s = X^i n i tT ahs(!) ' a ab s (i) an d n i is neutrino absorption cross section and num- 
ber density of target p article. Usual ly au thors define r < | or 1 as the criterion for neutrino 
transparent dCheng et al.l 120091 IJankal 1200 lb . Following we investigate chemical equilibrium condition 
for two different cases respectively. 



2.1 Case 1. Neutrinos are Transparent 

When the npe^ 1 gas is in equilibrium and transparent to neutrino and antineutrino (//„ = ^ = 0, that 
is, we hav e taken their number densities to be zero), the beta reactions are the most important physical 
processes dYuanl [20051) . The steady equilibrium state is achieved via the following beta reactions, 

e~ + p -> n + i/ e , (2) 
e + +n -> p + D e , (3) 
n ->• p + e~ + D e . (4) 

Reactions(f2]i-(|4li denote the electron capture(EC), positron capture(PC) and beta decay(BD) respec- 
tively. Since the system is transparent to neutrino and antineutrino, neutrinos and antineutrinos produced 
by reaction©-© can escape freely at once, inducing lots of energy loss, so the reverse reactions, neu- 
trino capture and antineutrino capture, are negligible. If the system is in equilibrium sate, composition 
is fixed and electron fraction Y e keeps as constant. EC decreases the Y e , while PC and BD increase Y e , 
then a general steady equilibrium condition is 

A e -p = A e + n + A„, (5) 

where A s are the reaction rates, subscript symbols denote reaction particles (the same in following 
section). Other reactions such as7 + 70e _ + e + o^ + P also exist, but they do not influence the 
electron fraction directly. These beta reaction rates can be obtained in the previous studies, we list them 
as below(We here employ the natural system of units with m e = ft = c = 1. In normal units, they 

would be multiplied by ( " ( l ^f c ) dYuanl 120051: bmganke & Martfnez-Pinedoll2000l) . 



A 
A 

2^ 



dE e E ePe (E e - Q) 2 F(Z, E e )f e , (6) 
dE e E ePe (E e + Q) 2 F(-Z,E e )f e+ , (7) 

A„ ~ A n „ J W dE e E ePe {Q - E e ) 2 F(Z + 1, E e )(l - / e ), (8) 

Considering charge neutrality, Y e = Y p , and the conservation of the baryon number, Y n +Y P = 1, so 
n p and n n in Eq.s (6)-(8) are equal to pY e NA and p(l — Y e )Na, respectively. F(±Z,E e ) is Fermi 
function, which corrects the phase space integral for the Coulomb distortion of the electron or positron 
wave function near the nucleus. It can be approximated by 

F(±Z, E e ) « 2(1 + s) (2p e i?) 2 (- 1 )e^| J ^±M| 2 , (9) 

here Z is the nuclear charge of the parent nucleus, Z = 1/0 for proton/neutron, s = (1 — a 2 Z 2 ) 1 / 2 , a is 
the fine structure constant, R is the nucleus radius, ?/ = ±aZE e /p e , T(x) is the Gamma function. We do 
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Table 1 steady state chemical equilibrium condition as neutrino is transparent for T — 10 9 K. fi p fi e 
and fj, n are chemical potentials without rest mass. 



Y e 


P 


p 


A e +n 


An 


I 

Pe 


Pn 


Pv 


C 




gem -3 


cm -3 s" 1 


cm" 3 s -1 


cm -3 s" 1 


MeV 


MeV 


MeV 




0.10 


1.50E+08 


8.56E+26 


3.52E+19 


8.56E+26 


0.84 


-0.43 


-0.62 


1.09 


0.20 


6.83E+07 


5.21E+26 


2.22E+19 


5.21E+26 


0.80 


-0.51 


-0.63 


1.07 


0.33 


4.27E+07 


3.72E+26 


1.63E+19 


3.72E+26 


0.78 


-0.56 


-0.63 


1.06 


0.40 


3.02E+07 


2.79E+26 


1.26E+19 


2.79E+26 


0.76 


-0.60 


-0.64 


1.05 


0.50 


2.29E+07 


2.14E+26 


9.98E+18 


2.14E+26 


0.74 


-0.64 


-0.64 


1.03 



Table 2 steady state chemical equilibrium condition as neutrino is transparent for T — 5 x 10 9 K. The 
notes are the same as those in Table 1. 



Y e 


P 


^e~~ p 


A e +n 


An 


i 

Pe 


/ 

Pn 


P P 


C 




gem" 3 


cm" 3 s" 1 


cm" 3 s" 1 


cm" 3 s" 1 


MeV 


MeV 


MeV 




0.10 


2.32E+08 


2.00E+29 


1.57E+29 


4.28E+28 


0.64 


-3.00 


-3.94 


1.95 


0.20 


8.38E+07 


9.51E+28 


7.73E+28 


1.79E+28 


0.45 


-3.49 


-4.08 


1.96 


0.30 


4.43E+07 


5.71E+28 


4.75E+28 


9.61E+27 


0.33 


-3.82 


-4.18 


1.97 


0.40 


2.71E+07 


3.71E+28 


3.14E+28 


5.64E+27 


0.23 


-4.10 


-4.27 


1.98 


0.50 


1.78E+07 


2.47E+28 


2.13E+28 


3.38E+27 


0.14 


-4.35 


-4.35 


1.99 



not adopt any limiting form for Fermi function. Comparing to the derivation of the rates in Yuan(2005), 
we consider additionally the Coulomb screening of the nuclei. f e and f e + are the Fermi-Dirac functions 
for electron and positron. f e = [1 + cxp ( Ee kT ' J " s f e + = [1 + cxp ( " Ee J r A ' e where k is the 
Boltzmann's constant, electron chemical potential p e can be calculated as follows(energy is in unit of 
m e c 2 and momentum in unit of m e c), 

Q_ /-0O 

P N A Y e = ^J o (fe- fe+)p 2 d P . (10) 

where A e = is the electron Compton wavelength. Note that the calculation method of chemical 
potential of electron (including the chemical potentials of proton and neutron in Eqs.(l 1)-(12)) also 
differs from the method in Yuan (2005). For one system with the given temperature T and densities p, 
electron fraction Y e can be determined by iteration technique of Eq.©. 

Figure 1 shows the T, Y e and p that satisfy the equilibrium condition. It can be found that the Y e 
decrease with the densities. As p > 10 11 gcm~ 3 , Y P tends to zero, e specially for the lower temperatures. 
This consists with the results in Fig. 5 of Ref. dReddv et al.lll998t) . At the high density, the (3 decay is 
almost forbidden and the positron capture rate is smaller than that of electron. In order to sustain the 
equilibrium, electron number density n e must be very low, which causes the Y e de crease obviously. Note 
that i t quite differs from the direct Urea process for strong degenerate baryons dShapiro & Teukolskyl 
fl983b . in which n p /n n > 1/8. The baryons here are nondegenerate since their chemical potentials 
(minus their rest mass) are very low, even minus. After p, T and Y e are found, chemical potentials 
/!„, p,p can be calculated as below(energy is in unit of m e c 2 and momentum in unit of m e c), 

P N A Y p = ^ 1^ ^ [1 + exp( ^_-_M£ )r i^ 

P N A (1-Y e ) = |^%2 [1+cxp( £n^ )r l dpi (12) 

where the conservation of the baryon number and the charge density are also included. 
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In order to describe the numerical relationship of fi e , fi n and fi p , we define a factor C: fx n — 
Hp + Cfi e . Table 1 and Table 2 are the results at T = 10 9 K and 5 x 10 9 K respectively. It can be 
seen from Table 1 that A e - p « A„ >> A e +„ , i.e., positron capture rate at this time can be ignored. 
C « 1 means that jit„ = fi p + \i e is valid. While from Table 2 one can find A e - p sw A e + n >> A„, i.e., 
beta decay becomes neglectable because lots of positrons take part in the reactions at high temperature. 
Correspondingly, the equilibrium condition becomes to fj, n = \i p + Cfx e , C ~ 2. It is quite different 
to the well known result M n = Mp + Me - This result was first observed by Yuan (2005), and a detailed 
explanation can be found in (lYuanll2005l) . Here we give a simple explanation that, A e - p oc n e n p oc f e f P , 
A e +„ oc n e +n n oc /„/ e +, so A e - p - A e+n oc / e / p - f n f e + = / p / e +(j^ - Considering / e « 



exp(- 



kT 

is valid as A 



),fe 



- w cxp ( e fe y° ), / P « exp ( p kl r p ) and /„ « cxp ( we find ^„ = ^ p + 2/i e 

= A e +„. For a more universal case, none of A e - p , A e + n and A„ can be ignored, the 
coefficient C will vary with the physical conditions. Figure 2 shows the coefficient C at different T and 
Y e . It can be found that C mainly depends on temperature T. When T < 10 9 K, C»l; when T from 
10 9 K increases to 5 x 10 9 K, C increases significantly from 1 to 2; when T > 5 x 10 9 K, C « 2. But 
when T > 3 x 10 10 K and Yg > 0.4, C is obviously larger than 2. The reason is that the fiducial analysis 
in reference (lYuanll2005l) ignoring the Fermi function. If we set Fermi functions are equal to 1, C rj 2 
is still valid. For the convenience to practical application, we give an analytic fitting formula that can 
facilitate application, 

C = 2-[l+cx P (^-^)]-\ (13) 

where A = [2.8643, 2.9249, 2.9785, 2.9902, 3.0094], B = [0.79138, 0.72181, 0.66331, 0.61813, 0.57999] 
corresponding to Y e = [0.1, 0.2, 0.3, 0.4, 0.5]. T 9 is the temperature in unites of 10 9 K (T 9 E [1 - 6]). 
The accuracy of the fitting is generally better than 1%. 



CD 
> 




11 14 



logp(g cm ) 



logT(K) 

Fig. 1 electron fraction Y e as a function of T and p for equilibrium state npe^ 1 gas. 
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Fig. 2 the coefficient C of chemical potential equilibrium condition fi n = fi p + C/j, e as a 
function of T and Y e . 



As an example, we introduce the application to electron fraction of GRB accretion disk. GRB is 
one of the most violent events in the universe, but its explosion mechanism is still not clear. Many au- 
thors support the view that GRB originates from the accretion disk of stellar mass black hole. Various 
accretion rates (from O.O1M s _1 to 10M Q s _1 ) bring quite significant difference to the disk structure 
and composition. As the temperature of the accretion disk is generally larger than 10 10 K, all nuclei are 
dissociated to the free nucleons, so npe gas can describe the composition well. For lower accretion 
rates (M < 0.1 Mq s" 1 ), the disk is t ransparent to neutrinos and an tineutrinos, and neutrino and an- 
tineutrino absorption are not important dSurman & McLaughlin[|2~004l) . Adopting the steady equilibrium 
condition, Y e of the disk model PWF99( lPopham et al. , 1999h (accretion rate M = 0.1M© s _1 , alpha 



viscosity a = 0.1, and black hole spin parameter a = 0.95) are obtained in Fig. 3. Dashed line and 
solid line are the result from the steady equilibrium condition and the full calculation by Surman et 
al. respectively. It shows that in the inner region of disk(from 20km to 120km), electron fraction from 
different methods are conform principally, which indicates the composition in the disk is in approximate 
equilibrium state, but our result is generally smaller than that of Surman et al. While in outer region of 
the disk, Y e deviates from equilibrium, and this deviation increases with the accretion disk radius. 

Surman and McLaughlin (2004) did not bother with specifying the radial profile of the temperature 
and the density of the accretion disk when calculating the electron fraction as a function of the radius 
for model introduced by Popham et al (1999). Here we rewrite the temperature and the density formulae 
of Popham et al. (1999)'s analytical model as 

T = 1.3 x 10 n a°- 2 A/f °- 2 ir a3 K (14) 
p = 1.2 x 10 u a~ 1,3 Mf 17 i?" 2 55 Migcm" 3 (15) 

where M\ is the mass of the accreting black hole in Mq, and R is the radius in gravitational radius r g 
(r g = GMi/c 2 , which is equal to 1.4767km for Mi = IMq). Since the explicit formulae are given, 
we adopt the equilibrium condition of npe ± gas to obtain some representative values of Y e in Fig. 4 at 
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the radius larger than the inner edge (6 gravitational radius) of accretion disk. One can find from Fig. 4 
that Y e have a rapid increase with radius because both densities and temperatures decrease rapidly when 
radius increase and the variation of Y e is very sensitive to density and temperature as shown in Fig. 1 . 
For the different accretion rate, the accretion rate is larger, the Y e is larger. This means the distribution 
of Y e along the radius highly depends on the structure equations of the disk. 



0.45 



£ 0.35 




0.25 



Fig.3 Y e as a function of accrection disk radius for model M = 0.1, alpha viscosity a = 0.1, 
and black hole spin parameter a=0.95. The dashed line shows Y e from steady equilibrium 
condi tion, while the solid line is the full calculation by Surman et al. dSurman & McLaughlinl 
12004 



2.2 Case 2. Neutrinos are Opaque 

In the neutrino-opaque and antineutrino-opaque matter, neutrino and antineutrino will be absorbed by 
proton and neutron except the reactions ©-(IUi as following, 



v. 



n 
P 



P, 
- n . 



(16) 
(17) 
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Fig. 4 Y e as a function of accrection disk radius for the thin disk analytical model(a=0. 1 , a=0, 
Mi=3). Long-dashed line, dotted line and dot-dashed line show Y e as the accretion rate M 
=0.01, 0.05, 0.1, respectively. Vertical solid line denotes the inner boundary of the accretion 
disk (6 gravitational radius). 



Table 3 steady state chemical equilibrium condition as neutrino opaque for T = 5 x 10 10 K. The notes 
are the same as those in Table 1 . 



Y e 


P 

g cm -3 


p 

cm -3 s" 1 


cm -3 s _1 


cm -3 s" 1 


cm -3 s" 1 


cm -3 s _1 


Me 

MeV 


MeV 


MeV 


C 


0.10 


2.32E+11 


7.92E+37 


7.89E+37 


7.88E+36 


7.56E+36 


1.23E+31 


10.18 


-15.14 


-24.65 


1.01 


0.20 


6.17E+10 


2.03E+37 


2.01E+37 


4.17E+36 


4.01E+36 


5.87E+30 


6.69 


-21.40 


-27.38 


1.01 


0.30 


2.46E+10 


7.35E+36 


7.26E+36 


2.48E+36 


2.39E+36 


3.10E+30 


4.37 


-25.95 


-29.60 


1.01 


0.40 


1.05E+10 


2.75E+36 


2.71E+36 


1.40E+36 


1.35E+36 


1.52E+30 


2.48 


-30.29 


-32.03 


1.01 


0.50 


3.40E+09 


7.56E+35 


7.40E+35 


5.59E+35 


5.43E+35 


5.17E+29 


0.75 


-35.94 


-35.93 


1.02 


By using sections 


abs _ JL(T? 




C +Q) 2 - 


i] 1/2 (i-/ e ; 


and (Tpkp = 




-Q)[(Eu e - 





Q) 2 — 1] 1 / 2 (1 — / e +), we obtain their rates(the natural units system) 

K c n = 2^L n nJ {E v + Q)[{E Vc + Q) 2 -1}^ 2 F(Z + 1,E^ + Q){1 - f e )E 2 v J v JE v ^ (18) 



A PeP = ^n p I (Eg - Q)[{E Pe - QY - lf^F{-Z+l,E Ve - Q)(l - f e+ )E 2 J 9 Jv e , (19) 

where j Ut and f„ e are the Fermi-Dirac distribution function of neutrino and antineutrino. f„ e = [1 
exp ( "'^t 1 "' fv e = [1 + cxp ( "",7^"" ■ The number densities of neutrino and antineutrino 
are 



kT 



in 



p 2 dp- 



1 



in 



-h^r^l + cxp h*J""i + exp (2^) 



p dp- 



1 



(20) 
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Table 4 the evolution of initial electron fraction at different steady state chemical equilibrium con- 
ditions, t is the time post bounce, R v is the neutrinospheric radius, L n is the number luminosity for 
neutrino and antineutrino, < E Vc > and < Ey > are the average e nergy of neutrino and antineutrino 
respectively. All parameters above refer to Ref.l Arcones et al., 2008T) 



t 


Rv 


T 


L n 


< E Ve > 


< E Be > 


P 


1 (:' 


Y b 

1 e 


C 


s 


km 


MeV 


lO^s" 1 


MeV 


MeV 


gem 3 








2 


10.55 


6.34 


6.05 


20.71 


25.64 


5.50E+11 


0.113 


0.084 


1.39 


5 


9.82 


5.14 


3.55 


17.1 


22.6 


1.30E+12 


0.050 


0.039 


1.22 


7 


9.68 


4.73 


3.03 


15.9 


21.69 


1.40E+12 


0.042 


0.035 


1.15 


10 


9.59 


4.37 


3.06 


15.05 


21.86 


2.00E+12 


0.029 


0.028 


1.03 



When n Vp = n Pe , i.e. number density of neutrino is equal to that of antineutrino, /i„ e = = 0. In this 
case, the equilibrium condition Eq.© becomes 

A e -p ~ A„ e „ = A e +„ — Ap eP + A n . (21) 

One can find from Table[3]that even at T — 5 x 10 10 K, C is still approximate to 1. In another word, for a 
system with neutrino and antineutrino are opaque and their chemical potentials are zero, /i„ = /j, p + /i e 
is always effective no matter what temperature is, just as expected. 

3 EQUILIBRIUM CONDITION OF NPE ± GAS WITH EXTERNAL NEUTRINO FLUX 

As discussed in Section 2, we only consider that npe ± gas is isolated, but for many astrophysical en- 
vironments, the external strong neutrino and antineutrino fluxes can not be ignored. These processes 
involve some complex and difficult problems that concern both the neutrino transport and the interac- 
tions with nucleons. Here we discuss the neutrino-driven wind (NDW) from proto-neutron star (PNS) 
as a typical example. NDW is regarded as the major site for the r- process nuc l eosyn t hesis according to 



as a typical example. JNDW is regarded as the major site tor the r-process nucleosynthesis according to 
the o bservations of metal-poor-star in the recent years (see e.g.. lOianl 120081 20001 iMartmez-Pinedol 
l2008h . Since the NDW is firstly proposed by Duncan et al. in 1986 dDuncan et all Il986l) . many 



detailed analysis for this process have been done by many authors, including Newtonian and gen- 
eral relativity hy drodynamics and the other physical inputs, e.g. rotating, magnetic field, term i nation 



shock and so ondOian & WooslevL 1 1996t IThompsonl 120031: iMetzger et all 120071: iKuroda et all 120081: 
iThompson et all l200lt iFischer et all 120091). A ba sic sc enario of r-process nucleosynsis in the NDW 
can be simply described as(see Martfnez-Pinedoi 120081) : soon after the birth of PNS, lots of neutri- 



nos are emitted from the surface of PNS; because of the photodisintegration of shock wave, the main 
composition at the surface of PNS is proton, neutron, electron and positron (i.e. npe ± gas); in the cir- 
cumambience of PNS, the main reactions are the neutrino or antineutrino's absorption and emitting by 
nucleons (so called ' neutrino heat region'); in the further region electron fraction Y e keeps as a constant 
and a particles are combined; above this region, other particles, such as 12 C, 9 Be, are produced till the 
seed nuclei; abundant neutrinos are captured by seed nuclei in succession. The p revious researches show 
that the steady state is a good approximation to the NDW in th e first 20 seconds dThompson et al.L [200 lb 
IThompsonl l2003HOian & WooslevL fl996t IFischer et aill2009l) : 

Usually, neutron-to-seed ratio, electron fraction, entropy and expansion timescale are four essen- 
tial parameters for a successful r-element pattern. It is very difficult to fulfill all those conditions self- 
consistently. Electron fraction Y e is one of the most important parameters. Recent research by Wanajo 
et al. shows that the puzz le of the excess of r-element of A = 90 may be solved if Y e can increase 
l-2% (IWanaio et alll2009t) . The evolution of Y e is usually obtained by solving the dif ferential equation 
group which is related to the EoS, neutrino reaction rates and hydrodynamic frame dThompson et all 
120011) . Initial Y e at the origin of wind is an important boundary condition. Considering the neutrinos are 
emitted from the neutrino sphere, Y e at neutrino sphere can be regarded as the initial Y e of the wind. For 
a given model, the initial Y e can be determined by making the assumption that the matter in neutrino 
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sphere is in beta equilibrium dArcones et all 120081) . To co mpare the results with the p revious work of 
Arcones et al., we employs the same PNS model M15-ll-rl dArcones et al I l2008ll2007h . The model has 
a baryonic mass of 1 .4 M Q , obtained in a spherically symmetric simulation of the parameterized 15 M Q 
supernova explosion model. Detailed research shows that there are a few a particles will appear at the 
neutrino sphere, but number density of a particle is much smaller than that of proton and neutron, so 
it is reasonable to ignore the a particle effect on electron fraction, i.e., the matter is regarded as npe 
gas. Simultaneity, although lots of neutrino and antineutrino are emitted from PNS, their number densi- 
ties are equal, which means fi Ve = /j,p e = 0. Since the neutrino and antineutrino are transparent to the 
matter at neutrino sphere, neutrino produced by reactions (f2]i-([3]) can not interact with nucleons, but for 
the neutrino and antineutrino come from the core region of PNS, absorption reactions ( fToT ) and (T% are 
permitted. Their rates are 

Xv ' n = ~ Ye)NA ' (22) 

A PcP = ^^af e s pP Y e N A , (23) 

where L n<v and L n ^ a are the number luminosity of neutrino and antineutrino respectively, R v is the 
neutrinospheric radius. Considering too many physical factors (EOS , transport equation and so on) will 
influence the number luminosity and the neutrino energy, we simply assume the number luminosity and 
the energy of neutrino and antineutrino are the same as those in the wind. Firstly, we obtain the electron 
fraction by using a general equilibrium condition \ e - p — X VeTl = A e +„ — X PcP + X n . In other words, if 
the density and temperature are fixed for the equilibrium system, the electron fraction is unique. Then 
the coefficient C in the chemical potential equilibrium condition is determined (leftmost column in 
Table |4]i . The results for model M15-ll-rl are sh own in Table H Yf" i s the electron fraction for an 
extreme case (7 = 1, which is adopted in reference dArcones et al 1 120081) : F e b is the result in which the 



steady equilibrium condition is valid and the external neutrino flux is also considered. We can find F e 6 is 
universal smaller than Y" e a , which means the external neutrino flux strongly influences the composition 
of equilibrium system. Comparing Y^ with Y%, one can find that the improved equilibrium condition 
makes the electron fraction decrease significantly when the time is less than 5 seconds post bounce. 
After 5 seconds the electron fractions are similar to the case (7=1. Note that it is just a conclusion 
for the model M15-ll-rl. Due to the huge difference between the different models, the results may be 
quite different for the other models. More detailed consideration will be done in our further work. Initial 
electron fraction is an important boundary condition to determine the electron fraction of the wind. Since 
r-process nucleosynthesis is strongly dependant on the electron fraction, the accurate electron fraction 
is useful for the final r-process nucleosynthesis. 

4 CONCLUSIONS 

In this work, we derive the chemical potential equilibrium conditions fi n = \i v + C\i e for npe^ 
gas at two cases ( with/without external neutrino flux). Especially in the neutrino-transparent matter, 
employing the fitting Eq.(|T3t for the transition from low temperature and high temperature is a 
more convenient method than the calculation of interaction rates as usual. Since chemical potentials 
are dependant on three parameters: density, electron fraction and temperature, any one of those 
three parameters can be determined if the other two parameters are given. Although the variation 
of factor C is complicated as the external neutrino flux cannot be ignored, one can obtain the ex- 
tremum of those parameters assuming the C = 1 or 2. Furthermore, our results can be regarded as 
the reference value for non-equilibrium sates. Considering the simplicity and the far-ranging astro- 
physical environment, the results in this paper is expected to be used widely in the further relative works. 
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